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1.0  Phase  II  Accomplishments 


During  Phase  II  our  team  has  tackled  the  main  challenges  related  to  the  demonstration  of  low  noise  and 
high-frequency  miniaturized  aluminum  nitride  (AIN)-based  frequency  sources.  Given  the  Phase  I 
demonstration  by  our  team,  we  have  especially  focused  on  understanding  the  noise  sources  in  our 
oscillators  and  identifying  new  methods  to  evade  phase  noise.  Briefly,  we  summarize  in  bullet  form  the 
main  accomplishments  and  contributions  made  in  understanding  and  engineering  miniaturized 
frequency  sources.  The  following  sections  of  this  report  provide  more  in  depth  details  relative  to  each 
accomplishment. 

1)  We  demonstrated  the  ability  to  use  parametric  bifurcation  to  filter  and  shape  the  noise  of  a  227- 
megahertz  (MHz)  AIN  oscillator.  The  resulting  phase  noise  (PN)  at  1  kHz  offset  scaled  to  a  1- 
gigahertz  (GHz)  carrier  is  -  111  decibel/hertz  (dBc/Hz).  The  obtained  PN  value  exceeds  the  Phase  II 
milestone.  A  general  theory  describing  this  phenomenon  was  also  derived. 

2)  We  identified  major  source  of  damping  in  high-frequency  AIN  contour-mode  resonators.  Anchor 
losses  are  dominant  for  frequencies  in  the  range  of  200-300  MHz,  whereas  interfacial  dissipation 
(electode/AIN  interface)  dominates  above  750  MHz.  Segmented  electrodes  were  introduced  to 
reduce  damping  without  affecting  the  resonator  transduction  efficiency. 

3)  We  demonstrated  the  ability  to  use  duffing  non-linearities  in  AIN  resonators  to  cancel  amplifier 
phase  noise  below  the  thermal  time  constant  of  the  resonator.  Practical  demonstration  is  described 
by  theoretical  model.  This  demonstration  can  be  broadly  applied  to  other  non-linear  resonators. 

4)  We  identified  frequency  flicker  noise  in  the  AIN  resonator  as  the  main  PN  limiting  mechanism  when 
a  low  noise  amplifier  is  used.  A  non-linear  model  describing  the  resonator  noise  was  generated  and 
imported  into  a  standard  circuit  simulator. 

5)  Temperature  compensation  via  resonator  ovenization  was  experimentally  demonstrated.  Two 
methods  of  using  ovenization  were  implemented  and  demonstrated.  A  final  temperature  stability 
of  <2  ppm  from  -  40  to  +  85  °C  was  demonstrated  for  a  microcontroller-based  ovenized  oscillator. 
Fundamental  limits  associated  with  resonator  aging  and  frequency  variations  were  identified  and 
will  be  considered  in  driving  Phase  III  developments. 

6)  An  acceleration  sensitivity  <  40  ppb/G  was  recorded  for  vibration  frequencies  between  100  Hz  and 
2  kilohertz  (kHz).  Acceleration  sensitivity  <  1  ppb/G  was  recorded  for  the  higher  vibration 
frequencies.  We  have  identified  the  major  sources  of  acceleration  sensitivity  to  be  in  the  component 
mounting  and  wirebonds. 

7)  The  temperature  and  acceleration  insensitive  devices  were  demonstrated  with  components  that 
occupy  a  total  volume  <  31  mm3.  The  majority  of  the  volume  is  occupied  by  the  discrete 
components  used  for  the  resonator  ovenization  circuit.  Plans  have  been  made  to  use  integrated 
circuits  for  such  a  demonstration  in  Phase  III. 
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2.0  Parametric  Bifurcation  Enables  Phase  Noise  Filtering 


We  excite  parametric  bifurcation  in  a  parametric  divider  formed  by  a  non-linear  varactor  and  two 
parallel  tanks  set  a  co0  and  2u>0,  u)0  being  the  oscillator  carrier  frequency  (Fig.  1).  The  parametric  divider 
is  placed  at  the  output  of  the  AIN  resonator  based  oscillator.  The  oscillator  output  frequency  is  first 
doubled  and  its  magnitude  amplified  prior  to  feeding  the  parametric  divider.  A  second  loop  taking  the 
output  of  the  parametric  divider  and  feeding  the  signal  back  to  the  amplifier  is  also  introduced  to  reduce 
the  amplifier  direct  current  (DC)-power  needed  to  activate  the  division.  A  parametric  divider  is  a  2-port 
network  that  transfers  energy  from  an  input  signal  at  co  to  an  output  signal  at  co/2.  Above  a  certain  input 
power  a  system  pole  becomes  unstable  and  a  flip-type  bifurcation  occurs.  The  unstable  pole  causes  the 
rising  of  a  sub-harmonic  oscillation.  This  activation  power  is  mostly  set  by  the  varactor  DC-bias,  Vdc,  and 
by  the  input  power,  Pin(2co0). 


$220  MH;  ■  Frequency  _  1,  . 

_  7.  i—*  Combiner 

Doubler 


Figure  1:  Parametric  Filtering  Topology  used  to  reduce  Oscillator  PN 

The  output  at  oj0  of  the  feedback-loop  oscillator,  based  on  a  227  MHz  AIN  resonator,  is  filtered  around 
the  oscillation  frequency  and  sent  to  a  frequency  doubler  (1).  The  output  of  the  frequency  doubler  is 
amplified  (2),  filtered  (3)  and  sent  to  the  parametric  divider  input  (4).  The  varactor  DC-bias,  Vdc,  and  the 
divider  input  power,  Pin,  regulated  by  the  amplifier  bias  voltage,  Vb,  sets  the  divider  working  point.  Half  of 
the  divider  output  is  sent  to  the  output  load.  The  other  half  is  sent  back  to  the  input  of  the  amplifier  (2)  to 

reduce  the  amplifier  DC-power  consumption. 


Considering  the  parametric  circuit  (Fig.  1),  and  using  the  voltage  across  the  varactor,  v(t),  as  the  only 
state  variable,  it  is  possible  to  define  a  complex  admittance  Yt(uj0),  function  of  Pin(2u)0)  and  Vdc: 

Yt(oj o)  =  YT(Pin( 2<u0),  2oj0,  Vdc  )  (1) 
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Under  the  assumption  that  YT(a;0)  is  real  close  to  the  division  threshold,  and  no  white  noise  is  added  by 
the  divider  circuit,  the  spectrum  of  the  PN  of  the  sub-harmonic  oscillation,  Acp^Q),  where  O  is  the 
frequency  offset  is  given  by: 


|  Ay>T(Q)|2= 


.  dYT 
d<p 

A 

dcp 


f+Af.tf 

dro 


(2) 


Eq.  (2)  shows  that  for  small  Q  the  PN  of  the  divider  tends  to  the  linear  prediction  (c(5s(Q)/2).  For  larger  Q 
(>  10-100  Hz)  and  especially  for  a  large  derivative  of  the  admittance  with  respect  to  frequency,  the 
second  term  at  the  denominator  plays  a  key  role.  By  regulating  the  varactor  DC-bias  it  is  possible  to 
modify  the  divider  working  point  and  make  it  operate  close  to  the  bifurcation  region,  where  dYT/dto  is 
maximized.  Effectively,  this  phenomenon  enables  reduction  of  the  oscillator  PN  because  the  divider  is 
not  able  to  track  the  fast  noise  perturbations  of  the  micro-electrical  mechanical  system  (MEMS) 
oscillator.  In  the  case  of  a  227  MHz  AIN  oscillator,  this  technique  brought  >  26  dB  improvement  in  the  PN 
at  Q  =10  kHz  offset  and  >  20  dB  improvement  at  Q  =1  kHz  offset  (Fig.  2-3).  This  technique  yields  a  PN  of 
-  124  dBc/Hz  at  1  kHz  offset  for  a  227  MHz  carrier.  Referenced  to  a  1  GHz  carrier,  this  PN  corresponds 
to  -  111  dBc/Hz.  To  our  knowledge  this  is  the  lowest  PN  ever  reported  for  any  MEMS  based  oscillator. 
The  demonstrated  technique  is  of  general  value  and  could  be  applied  to  any  MEMS  oscillator.  In  the 
appendix  we  provide  a  generalized  description  of  noise  filtering  that  occurs  in  a  parametric  system 
starting  from  the  Mathieu's  equation.  This  description  solidifies  our  experimental  demonstration  from  a 
theoretical  standpoint  and  offer  further  guidance  in  the  design  of  a  higher  frequency  oscillator  in  phase 
III. 


>26  dB 


N 

T  -loo 

Q 

CD 

D 

■150 


12  3  4  5  6  7 
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-to* 
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Figure  2:  In  purple:  PN  for  the  Feedback  Loop  227  MHz  Oscillator  without  Parametric  Filtering.  In 
black:  PN  at  the  Output  of  the  Parametric  Filter  when  Vdc  is  2.07  V  and  Vb=16V. 

In  this  case  the  system  is  working  close  to  the  division  bifurcation. 
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Figure  3:  Squares:  PN  at  1  kHz  offset  for  different  Varactor  DC-biases  and  for  Two  Distinct  Amplifier 
Bias  Voltages,  Vb  (Fig.  1).  Circles:  PN  at  10  kHz  offset  for  different  Reverse  Varactor  DC-biases  and  for 

Two  Distinct  Amplifier  Bias  Voltages,  Vb  (Fig.  1). 

In  the  cose  in  which  Vb=16V,  the  sub-hormonic  oscillation  occurs  for  a  much  narrower  range  of  varactor 
DC-biases  (2.07  -  2.15  V).  When  no  parametric  filtering  is  used  the  phase  noise  is  - 104  dBc/Hz  at  1  kHz 
offset  and  -125  dBc/Hz  at  10  kHz  offset.  Instead ,  the  use  of  parametric  filtering  enables  20  dB 
improvement  in  the  phase  noise  at  0=1  kHz  and  26  dB  at  0=10  kHz. 
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3.0  Mitigation  of  Anchor  Losses  and  Interfacial  Dissipation  in  AIN  Contour- 

Mode  Resonators 


Although  we  do  not  target  improving  the  resonator  Q  to  reduce  the  oscillator  phase  noise  (an  approach 
that  we  will  prove  of  little  effectiveness  in  the  section  on  resonator  flicker  noise),  it  is  important  to 
understand  the  sources  of  damping  in  our  resonator  to  ensure  that  Q  >  1,000  are  available  at  high 
frequencies.  This  is  essential  in  justifying  the  use  of  an  acoustic  resonator  over  lumped  components  (LC) 
in  conjunction  of  parametric  filtering,  as  it  will  yield  an  overall  lower  phase  noise  source. 

To  this  purpose  we  have  systematically  studied  the  behavior  of  324  devices  having  a  fixed  AIN  plate  size 
(60x144  pm),  operating  at  3  different  frequencies  (220,  370  and  1050  MHz),  and  having  variable  anchor 
size  (function  of  the  acoustic  wavelength.  A).  Each  particular  condition  was  tested  in  9  separate  devices 
so  that  meaningful  statistics  could  be  collected.  The  experimental  results  are  shown  in  Fig.  4.  It  is 
evident  that  anchor  losses  (energy  escaping  through  the  resonator  supports)  play  a  dominant  role  at 
frequencies  in  the  range  of  200-300  MHz.  At  higher  frequencies,  the  resonator  damping  is  instead 
dominated  by  interfacial  dissipation  (due  to  a  stress  jump  occurring  at  the  interface  between  the  AIN 
film  and  the  electrodes).  To  mitigate  both  sources  of  damping  in  the  resonators  we  have  included 
etched  air  slots  at  the  anchor  location  (Fig.  5-left)  or  introduced  segmented  electrodes  to  reduce  the 
electrode  area  coverage  without  impacting  the  device  electromechanical  coupling  (Fig.  5-right).  These 
two  approaches  have  yielded  respectively  a  50  %  improvement  in  Qat  220  MHz  (etched  slots)  and  at 
750  MHz  (segmented  electrodes). 


/=220MHz 


/=370MHz  /r=l. 05GHz 


Figure  4:  (Top)  Basic  Resonator  Geometry  and  Parametric  Variations  for  different  Frequencies. 
(Bottom-left)  Table  showing  average  Q  and  Standard  Deviation  (9  samples  each)  for  Resonators  of 
different  Anchors  Width,  Wa,  and  Length,  La.  (Bottom-right)  Plot  showing  how  Q  varies  largely  as  a 

Function  of  Anchor  Size  at  lower  Frequencies. 
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CONVENTIONAL  DESIGN 


SEGMENTED  ELECTRODES 


Figure  5:  (left)  Example  of  different  Types  of  Slots  that  were  introduced  in  the  Resonator  Body  near 
the  Anchor  in  order  to  minimize  the  Amount  of  Energy  Escaping  through  the  Supports,  (right) 
Comparison  between  a  Conventional  Resonator  Design  and  a  segmented  Electrode  One. 

The  segmented  electrode  design  reduces  the  electrode  area  coverage  without  altering  the  resonator 

electromechanical  coupling. 


3 . 1  Duffing  Resonator  Cancels  Amplifier  Noise 


In  phase  I  we  have  shown  that  the  resonator  self-heating  is  the  leading  source  of  a  softening  nonlinear 
behavior  in  our  AIN  contour  mode  resonators.  We  illustrated  that  the  measured  resonator  admittance 
exhibits  a  softening  behavior,  but  that  the  amount  of  bending  in  the  resonance  depends  on  the  rate  of 
frequency  sweep  employed  in  the  measurement.  We  further  showed  that  self-heating  in  our  AIN 
resonators  can  be  described  with  a  Duffing  resonator  model.  However,  this  analysis  was  performed  in 
the  quasi-static  limit.  In  this  phase,  we  considered  the  complete  dynamic  model  of  a  self-heating  AIN 
resonator  and  analyzed  the  response  of  the  resonator  to  a  frequency  sweep  through  resonance  and  its 
relationship  to  the  sweep  rate.  In  addition,  we  illustrated  that  while  the  self-heating  nonlinearity 
appears  similar  to  the  Duffing  nonlinearity  in  the  quasi-static  case,  its  dynamics  are  different.  This  has 
important  implications,  especially  related  to  the  ability  of  these  resonators  to  evade  amplifier  noise  in 
oscillators  using  the  techniques  described  by  Yurke,  et  at. 

We  introduced  a  non-liner  resonator  model  (Fig.  6),  which  describes  the  first  order  dependence  of  the 
mechanical  (motional)  components  (Rm,  Lm,  Cm)  on  temperature.  The  resonator  is  described  by  the 
following  set  of  differential  equations: 
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Figure  6:  (Top)  Microscope  Image  of  a  1  GHz  AIN  Contour-mode  Resonator  (Bottom)  Modified 
Butterworth  Van  Dyke  (MBVD)  Circuit  Model  with  Nonlinear  Motional  Elements 


By  solving  this  set  of  equations,  we  were  able  to  accurately  predict  the  resonator  response  when  excited 
by  an  external  network  analyzer  for  different  signal  sweep  rates  (Fig.  7).  These  results  imply  that  we  can 
correctly  picture  the  resonator  dynamics  and  use  the  aforementioned  model  in  a  closed-loop  oscillator 
to  predict  phase  noise  shaping  in  the  presence  of  strong  non-linearities. 
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Figure  7:  Experimental  (Discrete  Points)  and  Fitted  (Continuous  Line)  Response 

From  top  to  bottom:  admittance  magnitude,  admittance  phase,  and  resonator  temperature)  of  a  1  GHz 
AIN  contour-mode  resonators  for  different  input  signal  sweep  rates.  The  input  power  was  kept  constant 
at  11  dBm.  It  is  evident  that  the  resonator  behavior  depends  on  the  sweep  rate,  a  clear  sign  that  its 
dynamics  are  regulated  by  the  device  thermal  time  constant.  The  overlaid  dashed  black  line  represents 

the  fixed  point  solution  of  the  Duffing  resonator. 

Resonator  non-linearities  tend  to  worsen  the  oscillator  PN  by  converting  amplitude  noise  to  frequency 
noise.  Nonetheless,  this  non-linear  phenomenon  can  also  be  exploited  for  evading  amplifier  phase 
fluctuations  as  previously  shown  by  Yurke,  et  al.  In  fact,  a  non-linear  resonator  exhibits  specific 
phase/frequency  relationship  for  which  the  oscillator  dynamics  can  be  substantially  modified.  Intuitively, 
the  resonator  time  constant  can  be  made  infinite  at  special  operating  points  generated  by  the  resonator 
non-linearity.  When  the  resonator  is  locked  in  that  particular  state  in  a  closed-loop  oscillator,  then  the 
amplifier  phase  fluctuations  can  be  eliminated.  To  prove  this  theory  we  have  built  a  222  MHz  oscillator 
using  coaxial  electronic  components  (Fig.  8-top).  In  this  case,  the  AIN  contour-mode  resonator  (CMR)  is 
mounted  on  a  2-port  printed  circuit  board  (PCB)  (Fig.  8-bottom)  that  is  connected  in  series  with  the  rest 
of  the  components  in  the  circuit.  All  the  discrete  elements,  subminiature  version  A  (SMA)  connectors, 
and  cables  are  chosen  to  have  the  characteristic  impedance  equal  to  50  fi  for  matching  purposes.  An 
external  DC  power  supply  is  used  to  bias  the  radio  frequency  (RF)  amplifier  (Mini-circuits  ZKL-1R5+)  and 
the  output  is  taken  from  one  of  the  power  divider  (M/A-COM  T-1000)  terminals.  To  close  the  loop,  one 
terminal  is  connected  to  the  phase  shifter  (ATM  P1213).  In  this  particular  implementation,  the  selected 
amplifier  was  characterized  by  a  large  intrinsic  1/f2  noise.  When  driving  the  resonator  non-linearly  (by 
means  of  sufficient  gain)  and  tuning  the  loop  phase,  it  was  possible  to  identify  specific  working  points 
that  yielded  a  significant  phase  noise  reduction  (Fig.  9).  The  PN  reduction  was  limited  to  frequency 
offsets  below  the  resonator  thermal  time  constant.  The  obtained  results  can  be  fit  by  theoretical 
predictions.  It  is  important  to  note  that  PN  cancellation  occurs  only  for  noise  sources  external  to  the 
resonator.  When  the  noise  source  is  intrinsic  to  the  resonator  (see  next  section)  this  technique  is  not 
applicable.  Nonetheless,  the  practical  demonstration  is  of  general  interest  in  proving  the  use  of  non¬ 
linearity  for  shaping  phase  noise.  Our  team  has  ultimately  used  a  different  approach  based  on 
parametric  filtering  to  meet  the  stringent  phase  noise  requirements  dictated  by  the  DARPA  program.  For 
this  reason,  this  approach  will  be  considered  during  phase  III,  if  and  only  if,  ways  to  mitigate  the 
resonator  flicker  noise  will  be  identified. 
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Figure  8:  (Top)  Schematic  and  Actual  implementation  of  the  Oscillator  Setup  used  to  demonstrate  PN 
shaping  via  a  Duffing  Resonator.  (1)  AIN  CMR  mounted  on  a  PCB,  (2)  RF  Amplifier,  (3)  Power  Divider, 
and  (4)  Phase  Shifter.  (Bottom)  Photograph  of  222MHz  AIN  CMR  Wire  Bonded  on  a  2  port  FR4  Circuit 

Board. 
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Figure  9:  Measured  Phase  Noise  Response  of  the  Non-linear  222  MHz  AIN  Oscillator  for  different 

Phase  Conditions  (green  and  light  blue  curves) 

The  predicted  PN  response  is  overlapped  (blue  and  red  curves).  It  is  clear  that  the  oscillator  PN  is  reduced 
for  particular  amplitude  and  phase  conditions  exhibited  by  the  Duffing  resonator 


3.2  1/f  Resonator  Flicker  Noise 

In  order  to  identify  the  main  noise  sources  in  our  oscillators,  and  especially  explain  the  origin  (whether 
due  to  mechanical  or  electronic  components)  of  the  1/f3  slope  generally  seen  in  the  phase  noise  plots, 
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we  have  separately  measured  the  amplifier  baseband  noise  and  the  resonator  frequency  fluctuations. 
The  measurement  of  the  latter  was  complicated  by  the  lack  of  sensitive  equipment.  Ultimately,  we  were 
able  to  develop  a  method  to  measure  resonator  frequency  fluctuations  with  a  network  analyzer.  In  a 
simplified  way  the  method  is  based  on  observing  the  phase  fluctuations  in  the  reflection  coefficient  (or 
transmission  coefficient)  of  a  resonator  excited  by  a  fixed  frequency  signal.  As  long  as  the  resonator 
exhibits  a  noise  greater  than  the  signal  source,  it  is  possible  to  extract  the  resonator  frequency 
fluctuations.  We  have  measured  128  different  devices,  operating  at  different  frequencies,  having  varying 
AIN  film  thickness  and  electrode  materials,  and  fabricated  at  different  facilities  (University  of 
Pennsylvania,  Carnegie  Mellon  University,  Sandia  National  Laboratories  (SNL)  and  the  Institute  of 
Microelectronics  (IME)  in  Singapore).  Interestingly,  as  shown  in  Fig.  10,  there  is  no  dependence  of  noise 
on  device  Q,  volume  or  frequency  of  operation  (note  that  the  equivalent  closed-loop  noise  at  1  Hz  offset 
is  referenced  to  a  1  GHz  carrier).  It  appears  that  the  use  of  Pt  electrodes  reduces  the  resonator 
frequency  fluctuations.  Also  a  common  feature  for  the  devices  fabricated  at  SNL  and  IME  was  the  use  of 
an  oxide  film  to  isolate  the  resonator  from  the  high  resistivity  silicon  wafer.  Although  no  conclusive 
statement  can  be  made  at  this  stage  on  the  physical  source  of  frequency  fluctuations  in  the  resonators, 
we  can  confidently  state  that  the  oscillator  close-in  phase  noise  is  limited  by  the  acoustic  device  and  not 
the  electronics.  Fig.  11  is  a  clear  evidence  of  this  conclusion.  The  same  trend  was  also  seen  in  various 
other  oscillators. 
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Figure  10:  Plot  of  the  expected  Closed-loop  Phase  Noise  at  1  Hz  offset  referenced  to  1  GHz  Carrier  for 

128  different  Devices  vs.  Measured  Device  Q 

The  inner  color  of  each  symbol  corresponds  to  the  frequency  of  operation  indicated  on  the  scale  bar.  The 
outer  color  of  each  symbol  refers  to  the  exponent,  n,  of  the  phase  noise  (1/f). 
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Figure  11:  Comparison  of  measured  1  GHz  Oscillator  Phase  Noise  and  expected  Closed-loop  Oscillator 
Phase  Noise  based  on  Open  Loop  Measurement  of  Resonator  frequency  Fluctuations 

The  open  loop  measurement  was  performed  with  both  the  network  analyzer  and  the  Rhode  &  Schwarz 
FSUP.  The  two  instruments  give  the  same  expected  noise,  but  the  Rhode  &  Schwarz  FSUP  has  a  higher 

measurement  bandwidth  (up  to  1  kFiz). 

We  have  also  formalized  these  conclusions  by  generating  an  equivalent  circuit  model  in  advanced  design 
system  (ADS).  The  model  (Fig.  12-left)  describes  the  resonator  non-linearity,  its  frequency  flicker  noise, 
and  the  amplifier  (in  this  case  a  SiGe  BJT)  flicker  noise.  The  model  can  be  used  to  accurately  predict  the 
closed-loop  behavior  of  the  oscillator  phase  noise  (Fig.  12-right).  The  development  of  such  electronic 
model  is  a  key  step  forward  in  the  design  of  high  performance  MEMS  oscillators  and  will  have  a  broader 
impact  on  the  device  and  circuit  communities. 
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Figure  12:  (Left)  ADS  Oscillator  Circuit  Model  incorporating  AIN  CMR  Non-linearity,  Resonator 
frequency  Flicker  Noise,  and  Amplifier  Flicker  Noise  (Right)  Measured  PN  response  of  1  GHz  Oscillator 

compared  to  prediction  made  through  ADS  Model 

Note  that  no  fitting  parameters  are  used,  but  noise  sources  are  directly  measured  in  open  loop.  The 
overlap  between  simulation  and  experiments  confirm  the  validity  of  the  developed  model. 
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3.3  Temperature  Compensation 

In  Phase  II  we  used  the  same  temperature  compensation  approach  we  developed  in  Phase  I,  but  added 
the  ability  of  using  a  microntroller  to  fine  tune  the  oscillator  temperature  stability  and  meet  Phase  II 
metrics.  The  compensation  technique  is  based  on  integrating  a  resistive  heater  with  the  resonator  to 
either  ovenize  the  device  or  tune  its  center  frequency  in  response  to  temperature  changes.  As  shown  in 
Fig.  13,  we  looked  at  3  different  ways  of  integrating  the  heater  with  the  resonator:  (a)  under  the  body  of 
the  resonator,  (b)  on  the  top  surface  around  the  perimeter  of  the  resonator;  (c)  suspended  above  the 
resonator  body.  The  rationale  beyond  moving  the  heater  outside  the  resonant  body  was  to  limit  its 
impact  on  the  acoustic  performance  of  the  resonator.  This  is  done  at  the  expenses  of  increasing  the 
power  dissipated  to  elevate  the  resonator  temperature  because  a  greater  thermal  resistance  path  exists 
between  the  heater  and  the  device.  Ovenized  oscillators  were  built  using  all  3  different  heater 
implementations. 

Serpentine  on  the  Separate  serpentine  on 


bottom  electrode 


top  of  the  resonator 


Figure  13:  Three  different  ways  of  integrating  a  Serpentine  Heater  with  the  AIN  Resonator 

(Left  to  right)  heater  on  the  bottom  electrode,  heater  around  the  device  perimeter  of  the  resonator, 

heater  suspended  on  top  of  the  resonator 

As  in  phase  I,  the  use  of  a  Wheatstone  bridge  to  control  the  ovenized  oscillator  showed  a  temperature 
dependence  of  about  130-160  ppm  over  -  40  to  +  85  °C.  In  Phase  II,  attention  was  placed  in  selecting 
resistive  components  that  form  the  Wheatstone  bridge  with  very  low  temperature  coefficient  of 
resistance.  Nonetheless,  it  appears  that  the  temperature  stability  of  the  ovenized  oscillator  is  limited  by 
the  sustaining  amplifier.  A  major  improvement  from  phase  I  was  to  use  an  integrating  circuit  instead  of  a 
proportional  controller  to  stabilize  the  resonator  temperature  (Fig.  14).  Such  technique  enabled  us  to 
keep  a  more  stable  resonator  frequency  over  time. 
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Figure  14:  Comparison  of  Two  Ovenization  Techniques  to  achieve  Temperature  Compensation 

(Top)  proportional  circuit  used  in  Phase  I  based  on  a  Wheatstone  bridge ,  an  amplifier  and  a  switch  (BJT) 
to  control  the  resonator  temperature.  Although  an  overall  stability  of  about  130  ppm  was  attained,  the 
temperature  fluctuates  by  ~  ±20  ppm  over  time.  (Bottom)  integrating  circuit  implemented  in  Phase  II 
and  based  on  the  same  Wheatstone  bridge,  but  with  a  simple  amplifier  and  RC  network  in  the  feedback 
loop.  Similar  overall  temperature  stability  was  attained,  but  this  circuit  results  in  a  much  tighter  control 

of  frequency  over  time  with  fluctuations  <  ±  5  ppm. 

A  second  approach  that  uses  the  heater  as  a  mechanism  to  pull  the  resonator  frequency  enabled  us  to 
meet  the  program  requirements  (Fig.  15).  Frequency  stability  over  temperature  of  2  ppm  was  achieved. 
In  this  case,  the  resonator  center  frequency  is  programmed  to  be  at  a  particular  frequency  at  a  given 
temperature  by  storing  in  the  microcontroller  the  amount  of  current  that  needs  to  flow  in  the  resonator. 
This  method  is  agnostic  of  the  source  of  frequency  drift  (whether  the  resonator  or  the  circuit),  hence  it 
does  not  suffer  from  the  issues  of  the  Wheatstone  bridge  ovenized  oscillators  previously  presented. 

This  method,  though,  clearly  requires  a  calibration  process  that  measures  few  points  on  the  oscillator 
frequency-temperature  curve.  For  temperature  values  in  between  the  measured  ones  an  interpolation 
algorithm  is  used.  Given  the  use  of  unpackaged  devices  and  the  aging  of  the  resonator,  the 
microcontroller  needs  to  be  programmed  prior  to  any  testing. 


Time  domain  measurement  of  frequency  fluctuations  for  proportional  temperature  controller 

40r - 1 - 1 - » - 1 - 1 - 1 - 1 - 1 - 1 - 

If?  i|i  i»tij|ddiL|i'll tiJUbj  .  /)  1  ‘/I ' 

-1  *12349 


{»* 


13 

Approved  for  public  release;  distribution  unlimited. 


VddVdd 


Vheater 


Micro 

controller 


Temperature 

controller 


GaAs 

amplifier 


MEMS 


585,925,000 
585,924,000 

77  585,923,000 

z 

585,922,000 
§  585,921,000 
S’  585,920,000 
c5  585,919,000 
=  585,918,000 
6  585,917,000 
585,916,000 
585,915,000 

-30-20-10  0  10  20  30  40  50  60  70  80  90 

| _ Temperature  [°CJ _ 

Figure  15:  Microcontroller  based  Temperature  Compensation  Circuit 

An  external  digital  temperature  sensor  is  used  to  monitor  the  environment  temperature.  The 
microcontroller  is  first  calibrated  and  stores  the  amount  of  current  that  needs  to  be  flown  in  the  heater 
as  a  function  of  the  external  temperature.  A  practical  implementation  is  shown  in  the  middle  of  the 
figure.  A  GaAs  amplifier  is  used  to  drive  the  MEMS  resonator.  The  temperature  controller  is  mounted  on 
the  same  PCB  and  occupies  a  large  fraction  of  the  overall  volume.  A  temperature  stability  of  1.7  ppm  was 
obtained  over  the  temperature  range  of-  25  to  +  85  °C.  A  max  current  of  1.7  mA  was  used  with  a  heater 
supply  voltage  equal  to  12  V.  The  oscillator  was  not  tested  below  -  25  °C  to  avoid  electromigration  in  the 

Al  heater.  This  issue  is  addressed  by  using  a  Pt  heater. 


Ultimately,  we  believe  that  an  approach  based  on  the  resonator  ovenization  is  simpler  and  will  yield  a 
more  easily  field-deployable  component.  For  this  reason,  we  have  focused  on  finding  ways  to 
compensate  for  circuit  temperature  sensitivity.  By  properly  designing  the  bias  currents  that  feed  the 
sustaining  amplifier,  we  were  able  to  minimize  variations  of  the  circuit  operating  point  with  temperature 
(Fig.  16).  We  strongly  believe  that  this  approach  in  conjunction  with  the  use  of  integrated  circuits  will 
enable  us  to  meet  Phase  III  temperature  compensation  metrics  in  a  very  small  form  factor. 
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Figure  16:  Temperature  Insensitive  Oscillator  Circuit 

By  properly  setting  the  direction  and  magnitude  of  the  current  flowing  in  the  bias  circuit  over 
temperature  it  is  possible  to  make  the  amplifier  set  point  independent  of  temperature.  This  concept  was 
proven  experimentally  with  an  uncompensated  resonator.  The  circuit  temperature  sensitivity  went  from 

40  ppm  to  4  ppm. 


3.4  Acceleration  Sensitivity 

Limited  work  was  done  in  engineering  the  device  acceleration  sensitivity.  The  measurement  of  the 
acceleration  sensitivity  (Fig.  17)  clearly  highlights  its  dependence  on  factors  other  than  the  intrinsic 
resonator  out-of-plane  displacement.  The  oscillator  T  is  well  below  the  phase  II  metric  (<  30  ppb/G  for 
all  vibration  frequencies),  but  shows  a  linear  dependence  on  the  vibration  frequency,  which  is  likely 
indicative  of  poorly  mounted  components  or  wirebond  vibrations.  If  we  consider  the  vibration  frequency 
of  2  kHz,  we  can  see  how  our  existing  oscillator  already  meets  phase  III  specifications.  We  plan  to 
package  and  solder  our  components  with  a  more  permanent  fixture  so  as  to  extract  the  ultimate  device 
limit.  In  fact,  COMSOL  simulations  predict  that  the  device  itself  should  have  an  acceleration  sensitivity 
below  0.05  ppb/G. 
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Figure  17:  Acceleration  Sensitivity  Testing  performed  at  Vectron  on  a  Temperature  Compensated  Unit 

The  acceleration  sensitivity  is  <  30  ppb/G  for  all  vibration  frequencies.  It  is  clear  that  component 
mounting  plays  a  key  role  as  just  epoxying  some  of  the  components  reduces  spurious  peaks  in  the 

sensitivity  plot. 


3.5  Volume 

Phase  II  metrics  were  met  by  using  discrete  components.  Fig.  18  shows  the  overall  volume  taken  by  the 
components  mounted  on  the  PCB.  It  is  clear  that  the  greater  part  of  the  volume  is  consumed  by  the 
temperature  controller.  For  this  reason,  we  have  started  an  effort  on  using  integrated  circuit  (1C)  for  the 
implementation  of  the  temperature  controller.  This  controller  will  be  synthesized  and  tested  in  Phase  III. 


Oscillator  Circuit 


Component 

Area  [mm2] 

Thick,  [mm] 

Quantity 

Volume  [mm2] 

Transistor 

1.68 

0.7 

1 

1.18 

Resistor 

0.5 

0.5 

3 

0.75 

Capacitor 

0.5 

0.5 

4 

1 

Resonator 

0.25 

0.5 

1 

0.13 

Total  Volume  Oscillator  Circuit 

3.05 

Temperature  Controller 

Resistor 

0.5 

0.5 

3 

0.75 

Capacitor 

0.5 

0.5 

3 

0.75 

PIC12F1840 

9 

0.9 

1 

8.1 

Temp  sensor 

5.94 

0.9 

1 

5.35 

DAC 

5.94 

0.9 

1 

5.35 

Op.  Amp 

3.24 

0.8 

1 

2.59 

BJT 

5.94 

0.9 

1 

5.35 

Total  Volume  Temperature  Controller  28.2 

OVERALL  VOLUME  31.25 


Figure  18:  Volume  of  Discrete  Components  used  for  the  Synthesis  of  the  Ovenized  Oscillator 
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3.6  Phase  II  Metric  and  Comparison  with  PENNTAC  Team  Accomplishments 


Frequency 

[MHz] 

Phase  Noise  at  1  GHz 
for  1  kHz  offset 
[dBc/Hz] 

Acceleration 

Sensitivity 

[ppb/G] 

Temperature 

Coefficient 

[ppm] 

Volume 

[mm3] 

PHASE  II  METRIC 

800 

-110 

100 

10 

30 

PENNTAC  TEAM 

ACCOMPLISHMENTS 

>800 

-111 

<40 

2 

31 
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4.0  PHASE  III  PLANS 


Phase  III  goals  are  listed  below: 


Frequency 

[MHz] 

Phase  Noise  at  1 
GHz  for  1  kHz  offset 
[dBc/Hz] 

Acceleration 

Sensitivity 

[ppb/G] 

Temperature 

Coefficient 

[ppm] 

Volume 

[mm3] 

PHASE  II  METRIC 

1000 

-120 

10 

3 

1 

Briefly,  in  the  following  sections  we  outline  our  approach  and  risk  mitigation  for  each  metric.  In 
summary,  we  will  try  to  further  lower  the  oscillator  PN  by  using  the  parametric  filtering  technique 
identified  in  Phase  II.  Efforts  will  be  placed  in  miniaturizing  the  components  by  using  ICs  so  as  to  meet 
the  stringent  volume  requirements.  Low  noise  mechanical  and  electronic  methods  to  perform 
parametric  filtering  will  be  identified.  Simultaneously,  we  will  explore  techniques  to  mitigate  the 
resonator  intrinsic  noise  by  investigating  its  fundamental  physics  and  synthesizing  coupled  arrays. 
Finally,  we  will  identify  methods  to  further  reduce  its  temperature  sensitivity. 

4.1  Frequency 

This  is  a  low  risk  goal  for  our  team.  We  have  already  demonstrated  resonator  and  oscillators  operating 
up  to  1  GHz.  We  can  control  the  resonator  geometry  and  fabrication  to  produce  AIN  resonators 
operating  at  1  GHz.  We  are  also  familiar  with  the  required  electronics  to  design  the  electronic  oscillator 
at  this  frequency.  A  greater  challenge  is  posed  by  scaling  and  miniaturizing  the  electronics  required  to 
implement  the  parametric  filter  required  to  meet  the  Phase  III  metrics.  To  this  purpose,  we  have 
decided  to  take  an  approach  based  on  the  design  and  implementation  of  these  components  using  ICs. 

4.2  Phase  Noise 

This  is  the  highest  risk  metric  for  our  team.  Although  we  have  identified  a  pathway  that  simultaneously 
will  ensure  us  to  attain  low  phase  noise  and  miniaturization,  we  need  to  substantially  improve  the 
oscillator  phase  noise  by  another  10  dB.  We  plan  to  do  so  by:  1)  continuing  to  study  the  fundamental 
noise  sources  in  the  resonator  (which  limits  the  noise  in  the  sustaining  oscillator)  to  identify  the  physical 
origin  of  the  noise  and  methods  to  suppress  it;  2)  identifying  solid-state  or  mechanical  components 
capable  of  providing  parametric  noise  filtering  and  characterized  by  a  substantially  lower  intrinsic  noise 
than  off-the-shelf  components;  3)  using  coupled  resonant  systems  or  arrays  of  resonators  to  cancel  the 
resonator  intrinsic  flicker  noise  (Fig.  19). 
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Figure  19:  Examples  of  (top)  Dual  Mode  Resonator  and  (bottom)  Resonator  Array  that  will  be 
explored  in  Phase  III  to  reduce  the  Oscillator  Phase  Noise 

Addressing  the  issue  of  the  noise  source  in  the  resonator  is  also  essential  in  enabling  other  non-linear 
techniques  that  can  only  cancel  noise  sources  external  to  the  resonator.  We  believe  that  suppressing  the 
resonator  intrinsic  noise  could  enable  our  team  to  apply  other  noise  cancelling  techniques,  hence 
resulting  in  new  significant  gains  in  PN  shaping.  In  phase  II  we  have  developed  a  basic  theory  that  shows 
that  the  use  of  multiple  synchronized  references  can  be  exploited  to  cancel  the  intrinsic  noise  in 
resonators  (see  Q3-2013  report).  We  plan  to  implement  that  theory  in  Phase  III. 

4.3  Acceleration  Sensitivity 

This  is  a  low  risk  metric  for  our  team,  although  the  results  of  Phase  II  will  need  to  be  improved  in  order 
to  meet  the  final  goal.  Our  plan  is  to  identify  the  major  sources  of  low  frequency  vibrations  and  more 
appropriately  package  our  devices.  As  we  strongly  believe  that  our  oscillators  are  limited  by  external 
components,  we  are  confident  that  the  proposed  approaches  will  be  enough  to  ensure  attaining  an 
acceleration  sensitivity  <  1  ppb/G. 

4.4  Temperature  Coefficient 

This  metric  is  a  medium  risk  task  for  our  team.  We  believe  that  the  approach  used  in  Phase  II  can  be 
scaled  to  Phase  III.  Although  we  have  shown  that  temperature  stability  of  about  2  ppm  can  be  achieved, 
the  approach  is  limited  in  its  applicability  as  it  requires  a  very  stable  device  with  low  aging  and  drift. 
Resonator  aging  and  drift  directly  impacts  the  repeatability  of  the  measurements  and  limits  the 
maximum  stability  that  can  be  attained.  We  plan  to  study  the  resonator  aging  over  time  in  order  to 
bound  the  maximum  stability  that  can  be  achieved.  To  this  purpose  we  will  use  Vectron  packaging 
capabilities  and  study  our  resonator  behavior  over  periods  of  weeks.  We  will  also  investigate  the 
physical  reasons  of  aging  and  identify  methods  to  mitigate  it.  In  order  to  ensure  more  stable  ovenized 
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resonators  we  also  plan  to  introduce  compensation  layers  in  the  structure  of  the  resonator  so  as  to  start 
from  a  device  with  a  lower  intrinsic  temperature  coefficient  of  frequency  (TCF). 

Finally,  we  plan  to  implement  an  1C  version  of  the  temperature  controller  in  addition  to  a  temperature 
stable  gain  element  to  meet  the  volume  requirements  of  Phase  III.  To  this  effect,  efforts  were  already 
made  in  Phase  II. 

We  will  also  investigate  the  oscillator  start-up  time  according  to  MIL-PRF-55310E. 

4.5  Volume 

Phase  II  oscillator  volume  was  limited  by  the  size  of  the  temperature  controller.  In  order  to  meet  the 
phase  III  metric  we  plan  to  focus  on  the  miniaturization  of  the  temperature  controller.  To  this  purpose, 
we  plan  to  place  the  sustaining  oscillator  and  the  temperature  controller  on  the  same  chip  made  with 
integrated  circuits.  The  proposed  approach  will  result  in  an  overall  volume  of  1.12  mm3  (see  Fig.  20). 


Temperature  Compensated  Oscillator 


Component 

Area  [mm2] 

Thick,  [mm] 

Quantity 

Volume  [mm3] 

Resonator 

0.5 

0.25 

1 

0.125 

Oscillator  + 

TDC  and  DAC 

4 

0.25 

1 

1 

Total  Volume 

1.125 

Figure  20:  Expected  Oscillator  Volume  in  Phase  III 

The  miniaturization  of  the  temperature  controller  will  reduce  the  overall  volume  to  1.125  mm3. 
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5.0  APPENDIX 


The  Mathieu  Resonator  Model  of  Dynamic  Phase  Noise  Filtering 

A  simple  model  exhibiting  the  phenomenon  of  dynamic  phase  noise  filtering  is  the  Mathieu  resonator 
(our  parametric  divider). 


< 7  +  2  r  q  +  (cug  +  2o)0^(t)  +  45/72  cos(2/2t  +  20(t)  +  xpy)q 
—  V0  cos(/2t  +  0(t))  +  n(t), 


(Al) 


where  ^ (t)  is  an  intrinsic  fluctuation  in  the  resonant  frequency  of  the  resonator,  n(t)  is  an  additive 
white  noise  with  autocorrelation  (n(t)n(t'))  =  2 D8(t  —  t'),  22  is  the  mean  frequency  of  the  resonant 
drive  signal,  and  0(t)  is  the  phase  of  the  drive  signal  which  captures  the  phase  diffusion  and  frequency 
fluctuation  of  the  source.  We  assume  that  the  parametric  pump,  45/72  cos(222 1  +  2 0(t)  +  ip),  is 
generated  by  a  frequency  doubling  and  a  phase  shift  of  the  resonant  drive,  V0  cos(/2t  +  0(t))  (similarly 
to  our  experimental  setup).  Hence,  the  pump  signal  carries  the  phase,  2 0(t).  We  assume  that  the 
damping,  driving,  and  noise  n(t)  is  weak.  The  source  phase  noise,  0(t),  is  not  assumed  to  be  small,  but 
is  assumed  to  be  slow  such  that  0  is  small.  The  resonator  frequency  fluctuation,  ^ (0,  is  assumed  to  be 
weak  and  slow.  We  further  assume  that  the  mean  drive  frequency,  22,  is  close  to  the  natural  frequency 
such  that  H  =  co0  +  a  where  o,  the  frequency  detuning,  is  small.  Under  these  conditions  the  system  can 
be  analyzed  by  perturbation.  To  facilitate  this  we  first  introduce  the  coordinate  transformation, 

q  =  u  cos(/2 1  +  8)  +  v  sin(22t  +  0),  (A2) 

q  —  — (/2  +  8)u  sin(22t  +  0)  +  (/2  +  8)v  cos(22t  +  0),  (A3) 


where  u  and  v  are  the  new  state  variables.  Applying  equations  (A2)  and  (A3)  to  equation  (Al)  gives  a 
system  of  first  order  differential  equations  for  u  and  v  in  which  u  and  v  are  small  and  cyclo-stationary  in 
time.  This  makes  the  systems  amenable  to  the  method  of  averaging  [1],  [2].  The  perturbation  analysis 
results  in  the  set  of  slow  flow  equations. 


u«  —  r(l  +  5sini/>)u  +  r  a  +  8cosip  + 


8-J 

r 


V 


r  a  — 8 cosip  + 


0-<f 


v  +  n^t), 
Vn 


u-r(l-8  sini/Or  -  —  +  n2(t), 


(A4) 

(A5) 


where  %(/)  and  n2(t)  are  two  Gaussian  white  noises  with  (nj(t)n;(t'))  «  ^5i;-5(t  —  t').  We  assume 
that  all  fluctuations  in  this  system  have  zero  mean.  Thus,  the  underlying  average  flow  in  equations  (A4) 
and  (A5),  that  is  (u)  and  (i>),  has  a  fixed  point  at  ( u ,  v )  with. 


o  +  8  cos  ip  V0  _  1  +  5  simp  V0 

U  ~  52  —  cr2  —  1  2/72  ’  V  =  52  -  a2  -  1 2/72’ 
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(A6) 


The  eigenvalues  of  this  fixed  point  are  A+  =  — r(  1  ±  V52  —  cr2).  For  5  2  —  cr2  <  1  this  fixed  point  is 
stable.  At  |5|  =  VT  +  cr2  the  fixed  point  undergoes  a  bifurcation  and  becomes  unstable.  However,  we 
concern  ourselves  with  the  case  |5|  <  Vl  +  cr2.  Provided  the  fixed  point  is  stable  the  noise  terms  in 
equations  (A4)  and  (A5)  will  cause  the  system  to  fluctuate  about  this  fixed  point.  We  refer  to  this 
fluctuation  as  u  =  u  —  u  and  v  =  v  —  v.  This  fluctuation  is  assumed  to  be  small  and  so  we  disregard 
the  product  of  u  or  v  with  6  or  £  in  the  equations  for  u  and  v  to  arrive  at 


-1  -  8  sin  ip  a  +  8cosxp](u\  fu\  fn_p\,  (ni\ 
—a  +  Scosip  —  1  +  8 sin0j  \vJ  \vJ  ^  '  \n2J' 


(A7) 


It  is  interesting  to  examine  the  absolute  phase  of  the  resonator  response, 

v 

0  =  0  +  arctan— .  (A8) 

u 


For  weak  noise  the  fluctuation  in  the  resonator  phase  is. 


0  =  0-  0W0  + 


uv  —  uv 
u2  +  v2 ' 


(A9) 


It  is  straightforward  to  solve  equation  (A7)  and  calculate  the  power  spectrum  of  0.  The  expression, 
however,  is  rather  cumbersome.  We,  therefore  report  the  spectrum  only  for  the  particular  interesting 
case  in  which  the  resonator  is  driven  at  resonance  and  the  parametric  pump  is  phase  shifted  to 
attenuate  the  amplitude  of  oscillation:  0  =  7t/2  and  a  —  0.  In  addition  we  assume  that  0(t),  (( t ),  and 
n(t)are  uncorrelated.  In  this  case  we  have 

r2  l  4  D\ 

S4> (")  “  +  5)2  +  w2  [0-  +  8)2Se(a))  +  -  +  (i  -  5) 2 -^-J.  (A10) 

This  case  is  of  particular  interest  because  one  eigenvector  aligns  with  azimuthal  direction  in  the  ( u ,  v ) 
plane.  Indeed,  for  this  case  the  fixed  point  is  ( u ,  v )  =  (0,  —  y0/2/72(l  —  8 ))  and  the  eigenvectors  are 
e+  =  (1, 0)  and  e_  =  (0, 1)  for  8  <  0  or  e+  =  (0, 1)  and  e_  =  (1, 0)  for  8  >  0.  When  the  eigenvector 
is  orthogonal  to  ( u ,  v)it  is  aligned  to  the  azimuth.The  slow  eigenvector  is  in  the  azimuthal  direction 
when  8  <  0  and  the  fast  eigenvector  is  in  the  azimuthal  direction  when  5  >  0.  Note  that  8  >  0 
corresponds  to  parametric  amplification  and  additional  nonlinear  terms  can  become  important  as  the 
amplitude  of  oscillation  grows.  8  <  0  corresponds  to  parametric  attenuation  for  which  the  amplitude  of 
oscillation  is  suppressed  and  additional  nonlinear  terms  are  less  likely  to  become  important.  By 
manipulating  the  parametric  pump  coefficient,  8,  the  eigenvalue  along  the  azimuthal  direction  can  be 
made  slower  or  faster.  Now,  from  equation  (A10)  we  see  that  the  Mathieu  resonator  becomes 


22 

Approved  for  public  release;  distribution  unlimited. 


insensitive  to  the  phase  noise  of  the  source,  Sg,  as  5  -*  —  l.ln  this  limit  the  resonator  close-in  phase 
noise  will  be  dominated  by  S^/a)2,  which  is  identical  to  the  phase  noise  that  the  frequency  fluctuations, 
((t),  would  generate  if  the  resonator  described  by  equation  (Al)  were  to  be  used  as  the  frequency 
selective  element  in  a  closed-loop  oscillator.  Indeed,  if  we  seek  to  minimize  equation  (A10)  with  respect 
to  8  while  ignoring  the  white  noise,  /.e.let  D  ->  0,  we  find 


dS^ 

~d8 


2r2(l  +  8) 

CT2(l  +  5)2  +  oj2)2 


(sf(<u)-w2Sfl(<u)). 


(All) 


Equation  (All)  clearly  shows  that  when  the  frequency  fluctuation  of  the  source,  co2Se  —  Sg,  is  less  than 
the  frequency  fluctuation  of  the  Mathieu  resonator,  5^,  then  Sq is  reduced  by  increasing  8.  This  means 
approaching  the  bifurcation  point  at  8  ->  l.On  the  other  hand,  when  the  frequency  fluctuations  of  the 
source  are  greater  than  the  frequency  fluctuations  of  the  Mathieu  resonator,  then  S $  is  reduced  by 
decreasing  8,  thereby  approaching  the  bifurcation  point  at  8  -»  —1.  In  this  case,  .^approaches  the  limit 
set  by  the  Mathieu  resonator.  It  is  important  to  recall  that  this  model  is  valid  only  sufficiently  far  from 
the  bifurcation  point.  The  model  breaks  down  when  the  amplitude  of  oscillation  becomes  too  large  or 
the  system  is  placed  sufficiently  close  to  the  bifurcation  point  such  that  the  probability  for  the  noise  to 
push  the  system  out  of  the  domain  of  attraction  becomes  significant.  Both  considerations  require  the 
inclusion  of  nonlinear  terms  before  they  can  be  examined,  c.f.  [3],  [4], 

The  presented  theory  provides  a  generalized  description  of  PN  shaping  of  a  source  undergoing 
parametric  filtering.  The  agreement  between  theory  and  experiment  is  confirmed  by  how  the  PN  slope 
is  changing  when  filtered  by  a  parametric  resonator  as  shown  in  Fig.  A-l.  The  PN  changes  its  slope  from 
1/f3  (unfiltered  case)  to  1/f5  followed  by  regions  of  l/f3and  1/f2. 


Frequency,  Hz 

Figure  A-l:  Effect  of  Parametric  Filtering  on  the  Phase  Noise  of  an  External  Oscillator 

(As  described  by  the  theory  presented  in  this  section)  A  clear  change  in  the  slope  of  the  PN  is  observed 
when  filtered  by  a  parametric  resonator.  The  theoretical  prediction  is  in  line  with  our  experiments. 
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List  of  Acronyms,  Abbreviations,  and  Symbols 

ACRONYM  DESCRIPTION 

ADS  advanced  design  system 

AIN  aluminum  nitride 

CMR  contour-mode  resonator 

dBc  decibel 

DC  direct  current 

GHz  gigahertz 

Hz  hertz 

1C  integrated  circuit 

IME  Institute  of  Microelectronics 

kHz  kilohertz 

LC  lumped  component 

MEMS  micro-electrical  mechanical  system 

MHz  megahertz 

PENNTAC  Piezoelectric  Non-Linear  Nanomechanical  Temperature  and  Acceleration  Intensive  Clocks 

PCB  printed  circuit  board 

PN  phase  noise 

RF  radio  frequency 

SMA  subminiature  version  A 

SNL  Sandia  National  Laboratories 

TCF  temperature  coefficient  of  frequency 


28 

Approved  for  public  release;  distribution  unlimited. 


